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We consider how the Pade-Borel, Pade-Borel-Leroy, and conformal mapping sum- 
mation methods for asymptotic series can be used to calculate the dynamical critical 
exponent for homogeneous and disordered Ising-like systems. 

I. INTRODUCTION 

This paper shows in detail how the summation methods for asymptotic series can be 
used to calculate the critical exponent z determining the critical slowing down of the system 
relaxation time r ~ ~ |T — T^^^^ near the temperature Tc of the second-order phase 
transition ^ is the correlation length and u is the correlation length exponent). The most 
complicated and interesting direction in the theory of critical phenomena is the study of 
dynamical processes near a critical point. A significant achievement of the renormalization 
group (RG) approach in studying the static critical phenomena is the constructed concept 
of universality classes of the critical behavior of different systems characterized by similar 
critical properties. The universality of the dynamical critical phenomena in contrast to the 
static phenomena is considerably weaker, which is manifested in the existence of a wide 
range of models of the critical dynamics with different dynamical critical behavior P] and 
common critical properties in equilibrium. This difficulty is partially compensated because 
the dynamical critical characteristics of many of these models can be expressed in terms 
of the characteristics of their static critical behavior. This primarily concerns the models 
where the conservation laws hold for the order parameter or any other long-lived excitations, 
in particular, the energy density Q]. Because of the presence of preserved quantities, the 
purely hydrodynamic terms determining (and sometimes suppressing) the influence of the 
fluctuation processes on the dynamical (relaxation) behavior can appear in the stochastic 
equations of motion. The critical exponents describing the dynamical critical behavior can 
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be calculated in this case based on the static exponents with the same accuracy as the 
latter 

But there exist dynamical models such that the description of their critical behavior 
is not a trivial problem. They are primarily models with purely dissipative equations of 
motion. To study the critical behavior of such models, special methods must be developed, 
and dynamical quantities must be calculated independently of the static quantities. Among 
such models, the Ginzburg-Landau dynamical model (model A in the Halperin-Hohenberg 
classification scheme jl|) is most interesting; this model was first proposed by Landau and 
Khalatnikov to describe the anomalous sound attenuation in helium near a A-point. The 
dynamical critical behavior of other real systems such as the Ising-like magnets later came 
to be described in the framework of this model. 

The critical dynamics of model A was studied in P] using the RG methods based on 
the e-expansion. To describe this dynamics, one of us first developed a field theory ap- 
proach directly for the three- and two dimensional homogeneous systems successively in the 
three-loop and four-loop approximations and then also for structurally disordered 

three-dimensional Ising-like systems in the three-loop approximation [5] . But the dynamical 
critical characteristics were determined in these papers with an accuracy considerably less 
than the accuracy of the description of the static critical behavior of these systems in the 
six-loop approximation of the theory in j^Q,!^!^- This is primarily due to a fast increase in 
the volume of calculations already for the lowest-order perturbation series, especially in the 
description of the critical dynamics of structurally disordered systems. In this connection, 
the requirements for the accuracy of the summation methods for the asymptotic series in 
this theory used to calculate the value of the dynamical critical exponent z in the descrip- 
tion of the critical dynamics of model A in homogeneous and disordered systems become 
more rigid. Here, we first present the results obtained by summing the series in this theory 
using the Pade-Borel- Leroy (PBL) and conformal mapping (CM) summation methods to 
calculate the values of the exponent z. We compare the obtained results with experimental 
data and the results of a Monte Carlo computer simulation of the critical dynamics. 



II. MODEL 



The field theory RG method, which permits calculating the values of the critical expo- 
nents characterizing the asymptotic behavior of the thermodynamic and correlation func- 
tions near the critical temperature, is widely used to describe the anomalous properties of 
the thermodynamic characteristics of systems in the second-order phase transitions. Our 
model of the critical behavior of a homogeneous ferromagnetic system is the classical spin 
system thermodynamically equivalent to the 0(n)-symmetric Ginzburg-Landau model with 
the effective Hamiltonian 

n^lld'x{m' + ro^' + ^^'), (1) 

where d, is the dimension of the system, ip{x, t) is the n-componcnt order parameter (magne- 
tization) , ro ~ T — TcQ {Tco is the critical temperature in the mean- field approximation) , and 
g^o > is the vertex of the magnetization fiuctuation interaction. The dynamical behavior 
of a magnet in the relaxation regime near a critical point in the framework of model A is 
described by the Langevin-type kinetic equation for the order parameter, 

where Aq is the kinetic coefficient, C(^)^) is a Gaussian random force, and h{x,t) is an 
external magnetic field. It is well known that its solution in the form of correlation and 
response functions can be obtained using a generating functional of the form 



Q = y" exp ^-?^e//b, '^l + y" ('^^ + (p~h)d'^xdt^ , 

where an auxiliary field with the field source h and the action functional 



(3) 



Tie// = / d'xdt (Ao + ^^ (Ao^^ + ^) ) . (4) 

are introduced. The functions of the order parameter response to the field h are in this case 
determined as 



G{x,t) 



5 < ip{x, t) > 



5h{x, t) 



1 S^Q 



^<^{x,t)^iO,0)> . (5) 



^=0 ^5h{x,t)5h(0,0) 

Instead of the response function, it is more convenient to consider its vertex part r^^'^\k, u). 
The long-range long-lived fluctuations in the order parameter, which mainly determine the 



anomalies of the equilibrium and nonequilibrium characteristics of systems as T — T^, are 
taken into account in the framework of the RG methods developed in l|,llO|. The (A^+iV)th- 
order renormalized dynamical vertex functions r^^'^\k,uj), which are used in the theory 
and uniquely determine all the observable characteristics of the system, are given by the RG 
differential equation 



d d ^ d (N + N) 



rif = (6) 



with the renormalized charge g, the reduced temperature r, and the kinetic coefficient A. The 
asymptotic behavior of the thermodynamic and correlation functions, which have a power- 
law character as T ^ Tc, is determined by the zero g* of the function l3{g): (3{g*) = 0, and 
by the function values ■jrig*), l\{g*) and 70(5'*) specifying the static and dynamical critical 
exponents: z/ = (2 + 7r(5'*))^^, Tj = •y^^g*), z = 2 + 7a(5'*) (for example, if 7 = z/(2 — 77), 
then the magnetic susceptibility behaves as x ~ 1^ — Tc\~^). 

The functions l3{g), 7r(fi'), 1x{g) and 7(/,(5') contained in the differential RG equation can 
be calculated as series in ^f. If the space dimension d is close to four, then the coordinate 
of the fixed point g* of the function (3{g) takes small values. The methods of the pertur- 
bation theory in the coupling constant g A — d can be used in this case, and the critical 
exponents can be calculated. For real systems with d = 3,2 the series in g are asymptotic 
series, and they can be summed using special methods not based on the perturbation theory 
representations 0, Q]- 

In the description of the critical behavior of structurally disordered systems with frozen 
nonmagnetic admixture atoms or vacancies at the lattice nodes, which play the role of point 
defects of the structure, the additional term 

mv,V] = ^- I d'xV{x)^\ (7) 

where V{x) is the potential of the random field of defects, is introduced in the effective 
Hamiltonian of Ginzburg-Landau model (^, which leads to fiuctuations in the local critical 
temperature tq ~ T — Tco- The distribution of the structure defects over the volume of the 
system is assumed to be Gaussian (taking the deviations from the Gaussian distribution 
into account gives only corrections inessential in the critical domain) with the distribution 
function 



P[V] = Avexp 



(8) 



where Ay is the normahzation factor and Vq is a positive parameter proportional to the 
concentration of defects and to their squared potential. The correlation functions and the 
response functions for structurally disordered systems, obtained later using the generating 
functional, must be averaged over the random defect potential V{x). This averaging can 
be performed most effectively using the replica method, which is widely used for similar 
purposes in the study of weakly disordered systems (see, e.g., jlll). The essence of this 
method is that the following mathematical transformation is formally used to obtain the 
free energy of a disordered system in averaging over the defect distribution: 

~F/T = {{\nZ)) = \im-\n{{Z"^)) = lim ^ In Sp^^^^ exp(-^,,pi[{y,,}]/T), (9) 

where the double angle brackets denote the average over the probability distribution P[V] of 
different defect configurations, Z is the partition function of the original disordered system, 
and m replicas ("images") of the original field if - {^Pi} are introduced with the replica 
indices i = 1, . . . ,m. The replica Hamiltonian is then determined as 

exp(-i^S[{¥^.}]/T) = (|J[eM-H[{v^}]/T)jj , (10) 

and is translation invariant, in contrast to the original Hamiltonian. Applying this averaging 
procedure to dynamical generating functional 0, we can obtain the replica action functional 

+ 4t>o / d'^xdtdt'(pi{x,t)ipi{x,t)(fj{x,t')ipj{x,t'). (11) 

instead of the action functional for homogeneous model (j3)). According to the replica method, 
the properties of the original disordered system are obtained in the limit as m ^ 0. This 
limit annihilates the connected diagrams containing loops formed by the admixture vertex 

Vq. 

Taking the structural disorder of systems into account in describing the critical behavior 
thus leads to introducing an additional interaction vertex vq in the effective Hamiltonian. 
This vertex determines the effects of interaction of fluctuations of the nm-component order 
parameter in terms of the field of defects. The subsequent procedure for renormalizing the 
dynamical vertex functions is given by the RG differential equation 



d ^ d ^ d d ^ d (N + N) 
^dH + f^^d-g + ^^Yv - '^'d-r - ^^'dX 



rf'^^ = (12) 



with the renormahzation charges g and v. The functions f3g and (3^ in this case depend on 
g and v as on model parameters. The fixed points {g*,v*) of the RG transformations are 
determined by the zeros of the functions Pg and l3g{g*,v*) = 0, [3v{g*,v*) = 0, and the 
critical exponents are determined by the values of the functions 7^(5'*, f*), ■yip{g*,v*) and 
'j\{g*,v*) at the corresponding stable fixed points of the RG transformation. 

As in the case of homogeneous systems, the functions Pg, Py, 'jr, Jxid) and 70(5') can 
be calculated as series in g and v only if the system dimension d is close to four. For real 
systems with d = 3, the series in g and v are only asymptotically convergent. 



III. SUMMATION METHODS FOR ASYMPTOTIC SERIES 

To discuss the dynamical critical behavior of homogeneous systems and to calculate the 
exponent z, we only need the functions P{g) and 7a (5')- The explicit form of the first of them 
in the six-loop approximation for three-dimensional systems and in the four-loop approxi- 



mation for two-dimensional systems was obtained in 



We calculated the dynamical 



scahng 



unctions 7a in the four-loop approximation for two- and three-dimensional Ising 



3 ■ 



As a result, the functions P{g) and 7A(fi') can be represented as the series in the 



models 
charge g 

^ = -l+g- 0.71617%2 ^ 0.930767^^ - 1.582388^^ (13) 
9 

^^(g) = 0.027053^^ - 0.00418V + 0.022130/ 



for the two-dimensional model and 
^^^^ -l + g- 0.422497^2 + 0.351070/ - 0.376527^^ + 0.495548/ - 0.749689/, (14) 



9 

^^(g) = 0.008399^2 _ 0.000045/ + 0.020423/, 

for the three-dimensional model. Based on the method proposed by Lipatov [12], it was 
shown in that although the general term of the series for the function P{g) increases 
factorially, 

00 

/5(/ = J]c„/, Cn ^ c{-arn'n\[l + 0(l/n)], (15) 

n=0 

the series nevertheless satisfies the asymptotic convergence condition 

N 



P[g)-Y,cng' 



n=0 



<Ca^^'[iN + l)\]''\gf-^' (16) 



where C and a are some constants, in the wedge-shaped domain \arg g\ < ^o" on the complex 
plane. Matching asymptotic expansion ()15|) and the values of the first-order coefficients 
gives information about all the terms of the series and permits reconstructing the function 
I3{g) approximately, but this requires a special summation procedure for asymptotical! 
conver gent series. Special summation methods for such series were developed in 



151 Il6 l: the most effective of them are the Pade-Borel (PB), PEL, and CM methods. 
In ij the Lipatov technique for estimating higher-order terms in this theory was gen- 
eralized to the critical dynamics, and it was proved that the Lipatov asymptotic expansion 
also holds for the dynamical RG function 7a. Therefore, the above summation methods, 
which were successfully used to analyze the static RG functions, can also be used to analyze 
the dynamical RG functions. 

For the three-dimensional disordered Ising model, the obtained representations of the 
functions (3g, (3^ and 7a as series in the charges g and v have the forms 

(3g ^ 3 308 2 104 185 2 ,(«) i , 

— = -1 + 9 + -t; 9 QV + > Wlqv\ 



g " 2 729" 81 " 216 



i+i>3 



1^2 92 2 50 95 2^ V- . , 

— = — l + -9 + f 9 qv f + > h)'-qv\ (17) 

V 3 729 81 216 m ^ ' ^ > 

3 



i+i>3 



7a = ^ 7*,i^V = -^t; + 0.008400/ + 0.030862^i; + 0.053240^^ - 

2J=0 

-0.012642/ - 0.041167/?; - 0.152964^?;^ - 0.049995^;^ 

where the coefficients , 6 ■ ^ and 7, j- were calculated in and in the respective three- 
loop and six-loop approximations. 

As a preliminary illustration of the PB, PEL, and CM summation methods, we consider 
the exactly solvable problem of determining the ground state energy of the linear anharmonic 
oscillator with the Hamiltonian 

U = y'^x^^ gx\ (18) 



This permits detecting special features of the use of each of the summation methods and 
understanding how the accuracy of the used approximation depends on the number of terms 
present in the series. According to [3], the oscillator ground state energy Eo{g) can be 



represented as a perturbation series in the anharmonicity constant g whose first few terms 
have the form 

^ , , ^ „ . 3 21 2 333 o 30885 . 916731 , 65518401 « 

Eo(o) = > CnQ =l + -q 9 H H r 9 +•••, (19) 

Q\y) ny ^ ^gy ^54^ ^q24 ^ ^ 4095 ^ 32768 ^ ' ^ ^ 

n=0 

and exact values of Eq depending on g are given in . It was shown in that the general 
term of series ()19|1 increases factorially as in expression (fTH|) . 

The PB method is based on the idea that a series of form (fTIHl can be written as the 
integral 

/oo 
e''B{gt)dt, B{g) = Y,Bn9\ B^ = ^. (20) 

In this expression, B[g) is called the Borel image. According to asymptotic expansion 
(I15|l . the Borel image converges in a disk of radius 1/a. We then apply the Fade approxi- 
mation to B{g), which consists in using a rational function of the form 

[L/M] = (M > 1), (21) 

whose expansion in a Taylor series (in a neighborhood of the point g = 0) coincides with 
the expansion of the Borel image as closely as possible. A function of form (PT|) has L + 1 
coefficients in the numerator and M + 1 coefficients in the denominator. The entire set of 
coefficients is thus determined up to the common factor; to be definite, we usually set 60 = 1- 
As a result, we obtain the total set of L + M + 1 free parameters. This means that the 
coefficients of the expansion of the function [L/M] in a Taylor series must generally coincide 
with the corresponding coefficients of series ()19|) . If the number of terms N in the series is 
finite, then the approximants [L/M] must satisfy the condition L + M < N. 

A generalization of the PB method is the PBL method for which formulas (j^n|l can be 
written in the extended form as 

CO ^ 

/oo 

t'e-'B{gt)dt, B{g) = Bug"", B^ = (22) 

n=0 L[n + 0+ i) 

where b is an arbitrary parameter. In the simplest realization of the method, we consider 
the values 6 = and 6 = 1, for example, and then analyze the variations in the series 
approximations resulting from variations in the parameter b, which permits estimating the 
accuracy of approximations obtained by the PB and PBL methods. In our case where the 
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FIG. 1: Domain of the Borel image analyticity on the complex plane with the cut (—00, — 1/a) 
along the real axis (a) and its conformal mapping to the unit disk (b). 



PBL method is applied to the series determining the characteristics of the critical behavior of 
homogeneous and disordered systems, for a test series, we choose series fll9p corresponding 
to the exactly solvable problem of calculating the anharmonic oscillator energy with the 
asymptotic convergence of the series, which is similar to that of series (jlHj) and p4|l in the 
theory of critical phenomena. Therefore, applying the PBL method to series (fT!?|l and ()14|1 
in the theory of critical phenomena, we plan to use the values of the parameter h that lead 
to the best approximations of series (fT^ . 

We consider the main principles of applying the CM method to the Borel image (for 
6 = 0). According to asymptotic expansion (|T5|l . the Borel image is an analytic function in 
the complex plane of g with the cut from —1/a to —00 (Fig. ^). Because the integration in 
formulas (j^^ is performed along the entire real axis, we must analytically continue Bi^g) to 
arbitrary complex g beyond the convergence disk \g\ < 1/a, and this continuation is realized 
by the conformal map g = w{u) with 



u 



w[u] = —- 



u 



(23) 



a[l-uf' " (l + a^)i/2 + i' 
which takes the plane with a cut to the unit disk |m| < 1 (Fig. [T]d). The representation of 



1.40 
1.35 
1.30 
1.25 
1.20 
1.15 
1.10 
1.05 

0.0 0.2 0.4 0.6 0.8 1.0 9 

FIG. 2: Comparison of the exact values of the anharmonic oscillator ground state energy Eq 
(black squares) for different values of the anharmonicity constant g with the results obtained by 
the PB summation method for different [L/M] (solid lines), the PBL method with the parameter 
b = 2.221426 for the approximant [1/1] (white circles), and the CPB method with the constant 
a = 3 (stars). 

B{g) as a series in u thus leads to the convergent series for any g: 

oo oo 

B{g) = J2Bng''\,=^^^^ =^ B{u) = Y,Unu\ (24) 

ra=0 n=0 

Indeed, all possible singular points P,Q,R,... of the function B{g) lie along the cut, and 
their images P, Q, Q', R, R', ... lie on the boundary of the disk \u\ = 1. Hence, the second 
of series fl24|) converges for all |m| < 1. The relationship between the coefficients Un and B^ 
can be expressed as 

Uo = Bo, Un = J2B"^[-) ^n+Z-i {n>l). (25) 

These formulas solve the above problem as follows: the Borel image B{u{g)) converges 
for any singular points g = oo, g = —1/a and g = go with g^ G (— oo, — 1/a), and its 
coefficients Un are related are related to the original coefficients and the parameter b by 




TABLE I: Sums of squared deviations of the values of the oscillator energy Eq calculated using 
the PB, PBL, and CPB methods from the exact values of Eq (the gaps show that the Borel image 
has a pole when the corresponding approximant is used). 



[L/M] 


[1/1] [1/2] [2/1] [1/3] [2/2] [3/1] [1/4] [2/3] 




0.013923 - 0.002703 0.016687 0.000105 0.002535 


b 

o2 


2.221426 1.582184 1.466092 1.276639 1.131538 1.194466 2.016677 1.050884 

0. 441/04 O.O/yUzU 4.Uoo4o4 

n nnnni 1 n nnnnnR n nnnnn'? n nnnnn7 in^^ n nnnnm n nnnmn n nnnn/i/i 

U.UUUUll U.UUUUUD U.UUUUUo U.UUUUU/ o- iU U.UUUUUl U.UUUUIU U.UUUU44 

0.000010 0.000008 0.000007 


c2 

'-'CPB 


0.000175 0.041795 0.000836 0.000307 0.000192 - - 0.000011 


[L/M] 


[3/2] [4/1] [0/6] [1/5] [2/4] [3/3] [4/2] [5/1] 


5pB 


0.000001 0.004142 2.609852 0.031233 - - 0.000247 0.000833 


b 

'-'PBL 


1.005754 1.068052 $ 2.445345 0.970348 0.959599 0.964316 1.006032 

4.194370 2.003840 
3 • 10-^ 0.000004 0.000007 2 • 10"^ 2 • 10"^ 2 • 10"^ 0.000008 

0.000006 2 • 10-^ 


c2 

'-'CPB 


0.000064 - 3 • 10-'^ 5 • 10"^ 0.003932 0.000003 



linear transformation ()25j) (see formulas ()22|l ). Because we have |m| < 1 for any g, the Fade 
approximation can be effectively applied to the series in the variable u. This procedure is 
called the conformal Pade-Borel (CPB) method. 

In Fig. 121 we compare the approximations applied to the series for the anharmonic oscil- 
lator ground state energy with the exact values of Eq obtained using the PB method with 
different types of the Pade approximant [L/M] and the PBL and CPB methods (for a = 3 
[isl ) with the approximant [1/1]. Fig. |21 shows that the approximation error of the PB 
method increases strongly as g increases, although the situation becomes better as the num- 
ber N of the series terms taken into account increases and the diagonal and nearly diagonal 
approximants are used. But the PBL method for b = 2.221426 (in expression (j^Uj) ) and the 
CPB method even for the approximant [1/1] give results close to the exact values of Eq in 
the considered interval of the variable g. These results can be compared in accuracy only 



TABLE II: Values of the fixed point g* and the dynamical critical exponent z for the two- and 
three-dimensional Ising models obtained by the PB, PBL, and CPB methods. 





[L/M] 


[1/1] [1/2] [2/1] [1/3] [3/1] [1/4] [2/3] 


d=2 


PB g* 

z 

PBL g* 

z 

CPB g* 
z 


2.1971 - 1.8210 1.8271 1.7702 
2.1262 2.0868 2.0873 2.0820 
1.6978 1.7098 - 1.7174 1.6849 

2.0753 2.0764 2.0771 2.0742 
2.1448 - 1.9137 - 1.8193 
2.1270 2.1018 2.0922 


a=3 


PB g* 

z 

PBL g* 

Z 

CPB g* 

z 


1.5508 - 1.4314 1.4230 1.4175 1.4427 1.4753 
2.0202 2.0172 2.0170 2.0169 2.0175 2.0183 

1.3972 1.4019 1.4671 1.4056 1.4049 1.4073 

Z.UiOD Z.UiOO Z.UiOD Z.UiOD Z.UiOD 

1.5385 - 1.4552 - 1.4267 1.4253 - 
2.0452 2.0393 2.0374 2.0373 




[L/M] 


[3/2] [4/1] [1/5] [2/4] [3/3] [4/2] [5/1] 


d=3 


PB g* 
z 

PBL g* 

z 

CPB g* 
z 


1.4234 1.4292 - 1.4263 - 1.4200 1.3854 
2.0170 2.0172 2.0171 2.0169 2.0161 
1.4234 1.4355 1.4082 - 1.4221 1.4208 1.3848 
2.0170 2.0173 2.0167 2.0170 2.0170 2.0161 
1.4242 1.4245 - - - - 1.4147 
2.0372 2.0373 2.0366 



with the results obtained using the "best" diagonal approximant [2/2] in the PB method. 

In realizing the PBL method, it is more convenient to perform the integral Borel trans- 
formation as 

/oo 
dte-'B{gt''), B{g) = Y,Bn9\ = WTT^y (26) 

n=0 ^ + 

where the parameter b for all [L/M] is determined by the requirement that the mean squared 
deviation of approximations from the exact values of Eq be minimum on the entire interval 



of g. Using expression ()26|) allows approximating the exact value of Eq much better than 
using (|22I)- Moreover, this transformation, in contrast to transformation (j22I); allows avoiding 
the appearance of singularities in the integrand when the abovementioned variational method 
for choosing the values of h is used. 

Table HI presents values of the sum of squared deviations S"^ for the values of the oscillator 
energy Eq when different approximations of the exact values of Eq are used in the PB, PBL, 
and CPB methods. In Table HI we do not present the results obtained using approximants 
of the form [0/A^], because it was found that for even A^, using them leads to very large 
deviations of the series sum from the exact solution and the deviation S'^ does not have 
a minimum in h (only a small decrease in 5*^ is observed as h increases), while for odd 
N , singularities unremovable for any value of h appear in the Borel image corresponding 
to this approximant. For each of the approximants of the form — 1], we found two 
values of the parameter h characterized by close values of the deviation S"^. Comparing 
the given values of S"^ clearly shows that the accuracy of the PBL method is much higher 
than that of the PB method. The PB method permits obtaining results comparable in 
accuracy to the results obtained by the PBL method only starting from the fifth-order 
nearly diagonal approximants. This makes the PBL method preferable for analyzing the 
short series obtained in the description of the critical dynamics of homogeneous systems 
flT^ . (fT^ and disordered systems (fTTj) . We also note that the CPB method is second in 
accuracy among the methods considered and is as good in accuracy as the PBL method 
starting from the fifth-order approximants in N (except for the case of the approximant 
[4/2]). The PB method, although it is below the CPB method in accuracy for short series, 
permits obtaining comparable results only for the diagonal approximant [2/2] and the nearly 
diagonal approximants [3/2] and [4/2]. 

We now consider how the above methods can be used to calculate the dynamical critical 
exponent z for the homogeneous Ising model describing the critical behavior of systems with 
the dimensions d = 2 and d = 3. As previously noted, the exponent z is determined when 
the functions f3{g) and 7a(5') are used. The fixed point g* of the RG transformations is found 
from the condition l3{g*) = 0. In Table ITTl we present the values of the fixed point g* and the 
dynamical critical exponent z for the two- and three-dimensional Ising models, which were 
obtained by the PB, PBL (for the values of the parameter b related to the corresponding 

r □ 

approximants in the oscillator problem), and CPB methods, where, according to p,|2l, the 



constant a for the last method was taken to be 0.238659217 for d = 2 and 0.14777422 for 
(i = 3. In the PB and PBLmethods, the series for the function 2 + 7^ was summed using 
the fourth-order approximant [3/1], which is the best (summable) approximant for these 
methods, while in the CPB method, all the approximants other than [A^/0] turned out to 
be nonsummable, and only the simple CM method was therefore realized. 

The averaging procedure for the results obtained using the PB method with the approxi- 
mants of order > 4 taken into account allows obtaining the following values of the charge 
g* at the fixed point and of the exponent z: 

g* = 1.7987 ± 0.0201, z = 2.0847 ± 0.0019 {(1 = 2); 
g* = 1.4270 ± 0.0074, z = 2.0171 ± 0.0002 (d = 3), 

Using the PBL method gives 

g* = 1.7012 ± 0.0115, z = 2.0757 ± 0.0010 (d = 2); 
g* = 1.4125 ± 0.0046, z = 2.0168 ± 0.0001 {d = 3) 

Using the CPB method gives 

g* = 1.8193, z = 2.0922 {d = 2); 

g* = 1.4231 ± 0.0019, z = 2.0372 ± 0.0001 {d = 3). 

We note that all the obtained values of the charge g* at the fixed point for d = 3 agree 
well with the result g* = 1.416 ± 0.005 obtained by Le Guillou and Zinn- Justin 
Averaging the above values obtained by the different methods finally gives 

g* = 1.4209 ± 0.0035, z = 2.0237 ± 0.0055 {d = 3); 
g* = 1.7731 ± 0.0297, z = 2.0842 ± 0.0039 {d = 2). 

In the case where the critical behavior of disordered systems is described, the theory series 
for the functions Pg, /3„ and 7a (see (fTTj) ) are two-parameter series: 

F{v, 9) = Y. ^^^9' {i, J = 0, 1, 2, . . . ). (27) 

The methods described above cannot be applied to such series directly. They must be 
modified for the case of several variables. One that can be thus modified is called the 



TABLE III: Results obtained by the PB, PBL, and CPB methods in calculating the fixed points 
and the dynamical critical exponent z for the disordered Ising model 



[L/M] 


N=2 
[1/1] 


N=3 
[1/2] 


N=4 

[2/2] [3/1] 


N=5 

[1/4] [3/2] [4/1] 


N=6 

[4/2] [5/1] 


PB V* 

9* 
z 


-0.5868 
2.3900 
2.1509 


-0.6894 
2.2402 
2.1756 


-0.7383 -0.6991 
2.2743 2.2411 
2.1878 2.1780 


-0.6394 -0.7146 -0.7201 
2.2497 2.2640 2.2705 
2.1633 2.1819 2.1833 


-0.7125 
2.2583 
2.1814 


PBL V* 

9* 
z 


-0.6862 
2.3895 
2.1757 


-0.7001 
2.2499 
2.1787 


-0.7000 -0.7142 
2.2499 2.2459 
2.1787 2.1822 


-0.6900 -0.7142 -0.7089 
2.2601 2.2638 2.2646 

O i "71'' O r\ 1 OOO O 1 O A f\ 

2.1762 2.1823 2.1809 


-0.7141 
2.2600 
2.1822 


CPB i;* 

9* 
z 




-0.6800 
2.1654 
2.1716 


-0.6700 -0.6808 
2.1342 2.2562 
2.1691 2.1710 


-0.7302 -0.7239 
2.3079 2.2671 
2.1844 2.1831 


-0.7321 -0.7321 
2.2763 2.2799 
2.1854 2.1853 



A-method. The essence of this method consists in introducing a generahzed series 

oo 

F{v, g;X) = J2 ^r^iv, 9)X\ Cn{v, 9) = J2 C.,jV'9'8i+,,n, (28) 

n=0 i,j 

which converges to the original series (j?7|) for A = 1. We can assume that such a series 
depends on the single variable A and its coefficients are functions of u and v. The Borel 
transform of series (j27j) has the form In the A-method, the Borel image is reduced to a 
function depending only on the single variable A, the Pade approximant method is then 
applied to this function, and the integral is calculated. In the CM method applied to two- 
parameter series, we use the change 5 = v/u, which permits rewriting series (P7j) as 

oo 

F{v.9)\,=,. = Y.^n{5)v\ (29) 

n=0 

The Borel transform then becomes 

oo oo ~ 

F{v, 6) = I dte-'B{vt\ 6), B{v, 5) = ^ T{nb + if ^- (3°) 



The conformal map is determined by the formula 



uiv) = , = . (31) 

v/l + a(5)t; + l 



When the CM method is apphed to the series of the critical behavior theory for disordered 
systems, the problem becomes more complicated because the asymptotic behavior of the 
series is unknown in this case. In for the Ising model, it was proposed to take a{6) in 
the form 

a{6) =a{9/8 + 6) (5<0,5>4), 

where a = 0.14777422 for the three-dimensional Ising model. 

We now consider the application of these two methods to series (fT7j) . which characterize 
the nonequilibrium critical behavior of the disordered Ising model. In Table llllj we present 
the results obtained by these summation methods in calculating the coordinates of the fixed 
point, which determines the critical behavior of disordered systems, and the exponent z. In 
the PB and PBL methods, we summed the series for the function 2 + 7^ using the third-order 
approximant [2/1], which is the best for these methods, and in the CPB method, we used 
the approximant [1/2]. It follows from Table ITTTl that all these methods give close values for 
the vertices at the fixed points and for the exponent z. 

We now apply the averaging procedure to the obtained results using all the values of 
the approximants [L/M] except [1/1]. The averages of the fixed-point coordinates and the 
exponent z obtained using the PB method are 

t;* = -0.7019 ±0.0111, ^* = 2.2569 ±0.0048, ^ = 2.1788 ± 0.0027, 
Using the PBL method gives 

^;* = -0.7059 ± 0.0033, / = 2.2563 ± 0.0026, 2 = 2.1802 ± 0.0008 

Using the CPB method gives 

= -0.7070 ±0.0100, ^* = 2.2410 ±0.0227, z = 2.1786 ± 0.0026. 

The close values of the fixed-point coordinates and the exponent z obtained by different 
methods and practically coinciding within the accuracy confirm the reliability of the obtained 
results. Averaging these values obtained by different methods finally gives 

= -0.7049 ±0.0013, c/* = 2.2514 ±0.0042, 2 = 2.1792 ± 0.0004. 

The dynamical critical exponent z was calculated in Q using the Chisholm-Borel method 
and z = 2.1653 was obtained. This value of the exponent agrees rather well with the results 
obtained here. 



IV. CONCLUSION 



In this paper, we first successively applied the PB, PBL, and CM summation methods 
with subsequent use of the Fade approximation (the CPB method) to determine the values 
of the dynamical critical exponent z for homogeneous two- and three-dimensional systems 
described by the Ising model and for disordered Ising-like systems. Our analysis was based 
only on the now available short series of the theory corresponding to the results of the 
four-loop description of the critical dynamics of homogeneous systems and the three-loop 
description of disordered systems [J, 1^ . 

We compare the obtained results with the experimental data and the results of Monte 
Carlo computer simulations of the critical dynamics. In l2o(, the critical dynamics of the 
homogeneous Ising antiferromagnet FeF^ was investigated experimentally, and the value 
z = 2.1 ± 0.1 was obtained, which in the framework of a sufficiently wide interval of mea- 
surement error for z does not contradict our results of the field theory calculations using the 
summation methods. The Monte Carlo numerical investigations of the critical dynamics of 
the homogeneous three-dimensional Ising model for different data given in the literature gave 



the following values o: 



221, z = 2.04 ±0.03 



the dynamical critical exponent: z = 1.99 ± 0.02 12111 . z = 1.97 ±0.08 



2.04 ±0.01 



2^ and z = 2.032 ± 0.004 |25|. These values of 



the exponent z agree well with the results obtained here and confirm that the our predicted 
value z = 2.0237 ± 0.0055 is sufficiently reliable. 

As for the two-dimensional homogeneous Ising model, the numerical investigations of 
its critical dynamics, in contrast to similar investigations of the three-dimensional model, 
are characterized by a much wider range of the values obtained for the dynamical critical 
exponent: z = 2.14 ± 0.02 Q, z = 2.13 ± 0.03 ^ = 2-076 ± O.OOsQ, z = 2.24 ± 0.04 



2.24±0.07 ,3Ql, z = 2.16±0.04 



3l| and 2 = 2.1667±0.0005 |32j. This range of 2 is 



wide mainly because the critical behavior of the two-dimensional Ising model is characterized 
by both much greater amplitudes of the order parameter fluctuations in the critical domain 
(compared with the three-dimensional model) and more significant effects of the critical 
delay. The high-temperature expansion method used to describe the critical dynamics of 
the two-dimensional Ising model gives the value z = 2.125 js^. The values of the exponent 
z obtained for the two-dimensional Ising model in the cited papers are in the rather wide 
interval 2.07 < z < 2.31, while the values obtained here are on the lower boundary of this 



interval. This makes the resuhs obtained here more valuable as reference points in future 
studies. 

We now compare the values of the dynamical critical exponent for disordered systems 
calculated in this paper with the results of computer simulations of the critical dynamics 
of the three-dimensional disordered Ising model: z = 2.19 ± 0.07 for systems with spin 
concentration p = 0.95, z = 2.20 ± 0.08 for p = 0.8, z = 2.58 ± 0.09 for p = 0.6 and 
z = 2.65 ± 0.12 for p = 0.4 0; z = 2.16 ± 0.01 for p = 0.95, z = 2.232 ± 0.004 for p = 0.9, 
z = 2.38 ± 0.01 for p = 0.8 and z = 2.93 ± 0.03 for p = 0.6 Q- These results of critical 
dynamics simulations agree rather well with the results of field theory calculations using the 
summation methods only for weakly disordered systems with p > 0.8, while a noticeable 
difference between the results is observed for strongly disordered systems. We note that 
the results of the RG description of the critical behavior of disordered systems in this case 
hold only in the domain of weak disorder. To explain the dependence of the exponent z 
on the value of the structural disorder observed in computer simulations, the hypothesis of 
graduated universality was proposed in 22], according to which several types of different 
critical behavior can be observed in systems with spin concentrations greater than the spin 
percolation threshold depending on whether there exists only one spin percolation cluster in 
the system as in the case of weakly disordered systems or there is an admixture percolation 
cluster in addition to the spin percolation cluster as in the case of strongly disordered systems 
with transient regimes between the domains. 

Starting from the concept that the critical behavior of disordered systems is universal 
and that the asymptotic value (as L oo) of the exponent z is independent of the degree of 
disorder, the author of j3| obtained the asymptotic value z = 2.4 ± 0.1 using the effective 
values of the exponent listed above. But this value of the exponent is strongly inconsistent 
with the results obtained here. A numerical study of the critical dynamics of the three- 
dimensional Ising model with the spin concentration varying in a wide interval was analyzed 
in ^|. Assuming that the critical behavior of disordered systems is universal, the authors 
of 2^ obtained the asymptotic value of the exponent z = 2.62 ± 0.07 taking the effects of 
influence of the leading corrections to the scaling dependence for the dynamical susceptibility 



of the system into account. I n t 
uj = 0.50 ± 0.13, obtained in 



lis case, the value of the exponent of the scaling correction, 
is strongly inconsistent with the results of the field theory 



calculations of the static critical exponents, which were done using the summation methods 



in 3 



and gave uo = 0.25 ± 0.10. This value is also inconsistent with the results of numerical 



36| also done taking the effects 



investigations of the static critical behavior of the same model 
of influence of the leading corrections to the scaling dependence of thermodynamic quantities 
and the correlation functions with u = 0.37 ± 0.06 into account. In the approximations 



realized in the results for weakly disordered systems were characterized by the largest 
errors. We also note that the value of the dynamical critical exponent z obtained in 35| is 
even more strongly inconsistent with our results here. 

The critical dynamics of the disordered Ising antiferromagnet Feo.46-Z^^o.54-^2 was studied 
experimentally in [2^, where z = 1.7 ±0.2 was obtained. This result conflicts with both the 
results of the field theory calculations using the summation methods and the above results of 
computer simulations. This inconsistency can possibly be explained by a high concentration 
of nonmagnetic admixture atoms in the sample under study and by the presence of large- 
scale inhomogeneities in this sample, which significantly affects the characteristics of the 
nonequilibrium critical behavior. But in the other experimental paper concerning the critical 
dynamics of weakly diluted Ising magnet Feo.9Zno,iF2 [37], the value of the dynamical 
critical exponent z = 2.18 ± 0.10 was obtained as the result of high-precession measurement 
of the dynamical widening of Mossbauer lines using the Mossbauer spectroscopy method. 
This value agrees very well with the results of our calculations. Nevertheless, there is a strong 
need for additional investigations of the critical dynamics of disordered Ising-like systems, 
both experimentally and numerically. The results of calculations of the dynamical critical 
exponent for disordered systems obtained here will be reference points in such investigations. 
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